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ABSTRACT 

Context. Wide-field Spitzer surveys allow identification of thousands of potentially high-z submillimeter galaxies (SMGs) through 
their bright 24 yum emission and their mid-IR colors. 

Aims. We want to determine the average properties of such z ~2 Spitzer-selected SMGs by combining millimeter, radio, and infrared 
photometry for a representative IR-flux {A^^-st ~ 8 fim) limited sample of SMG candidates. 

Mefhods. A complete sample of 33 sources believed to be starbursts ("5.8 yum-peakers") was selected in the (0.5 deg-) J1046+56 field 
with selection criteria FjAfim > 400 yujy, the presence of a redshifted stellar emission peak at 5.8 yum, and r[,^^.^ > 23. The field, part of 
the SWIRE Lockman Hole field, benefits from very deep VLA/GMRT 20 cm, 50 cm, and 90 cm radio data (all 33 sources are detected 
at 50 cm), and deep 160;um and 70 ^m Spitzer data. The 33 sources, with photometric redshifts ~ 1.5 - 2.5, were observed at 1.2 mm 
with IRAM-30m/MAMBO to an rms ~0.7-0.8mJy in most cases. Their millimeter, radio, 7-band Spitzer, and near-IR properties 
were jointly analyzed. 

Results. The entire sample of 33 sources has an average 1.2mm flux density of 1.56 + 0.22 mJy and a median of 1.61 mjy, so the 
majority of the sources can be considered SMGs. Four sources have confirmed 4 cr detections, and nine were tentatively detected at 
the 3cr level. Because of its 24pm selection, our sample shows systematically lower -Fi.2mm/^24;im flux ratios than classical SMGs, 
probably because of enhanced PAH emission. A median FIR SED was built by stacking images at the positions of 21 sources in the 
region of deepest Spitzer coverage. Its parameters are T^^st = 37 + 8 K, Lfir = 2.5 x 10'' Lq, and SFR = 450 Mq yr"' . The FIR-radio 
correlation provides another estimate of LpiR for each source, with an average value of 4.1 x 10'^ Lq; however, this value may be 
overestimated because of some AGN contribution. Most of our targets are also luminous star-forming BzK galaxies which constitute 
a significant fraction of weak SMGs at 1.7 < z < 2.3. 

Conclusions. Spitzer 24 yum-selected starbursts and AGN-dominated ULIRGs can be reliably distinguished using IRAC-24 pm SEDs. 
Such "5.8 yum-peakers" with F24^ni > 400yuJy have LpiR > 10'^ Lq- They are thus z ~ 2 ULIRGs, and the majority may be considered 
SMGs. However, they have systematically lower 1.2 mm/24 ^m flux density ratios than classical SMGs, warmer dust, comparable or 
lower IR/mm luminosities, and higher stellar masses. About 2000 - 3000 "5.8 yum-peakers" may be easily identifiable within SWIRE 
catalogues over 49 deg- . 

Key words. Galaxies: high-redshift - Galaxies: starburst ~ Galaxies: active - Infrared: galaxies - Submillimeter ~ Radio continuvmi: 
galaxies 



1. Introduction 

Ultra-Luminous InfraRed Galaxies (ULIRGs, with 
^FiR <: lO'^Lo) are the most powerful class of star-forming 



* color figures and the Appendix are only available in the electronic 
form via http://www.edpsciences.org( 



galaxies. For 25 years, these prominent sources and their 
intense starbursts have been the tar get of ma ny comp rehen- 
sive studies, both locally (e.g., Sa nders '& Mirabel 19961 : 
iLonsdale et al. 2006t fV eilleux et al. 2009) and at high red - 
shift (e.g.. iBlain et aL 2004t iSolomon & Van den Bout 2005h. 
While local ULIRGs are relatively rare, submm/mm sur- 
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veys with large bolometer arrays such as JCMT/SCUBA 
[James Clerk Maxwell Telescope/Submillimetre Common 
User Bolometer Array (Holland et al. 1999)], APEX/LABOCA 
[Atacam a Pathfinder Experiment/Large Apex Bolometer 
Camera jSiringo etal. 2009 )1 or IRAM/MAMBO [Institut de 
Radioastronomie M illimetrique/Max-Planck Bolometer Array 
dKrevsa et al. 1998 1) 1 have shown that the comoving density of 
submillimetre galaxies (SMGs), which represent a significant 
class of high-redshift (z ~ 1-4) ULIRGs, is about a thousand 
times greater than th at of ULIRGs in the lo cal Universe (e.g., 
iLe Floc'h et al. 20051; IChapman et al. 2005h . They represent 
a major phase of star formation at early e pochs and are als o 
characterized by high stellar masses (e.g., iBorvs et al. 2005h . 
They are thus ideal candidate s to be the precursors of local mas- 
sive elliptical galaxi es (e.g., iBlain et aL 20021: iDve et al. 20081 
iLonsdale et al. 20091 hereafter Lo09, and references therein). 
Nearly all of the enormous UV energy produced by their mas- 
sive young stars is absorbed by interstellar dust and re-emitted 
at far-infrared wavelengths, with their far-infrared luminosity 
(Lfir) able to reach 10'^ Lq. However, despite the considerable 
efforts invested in mm/submm surveys, the total number of 
known SMGs remains limited to several hundred, and current 
observational capabilities are still somewhat marginal at many 
wavelengths. We thus still lack comprehensive studies of SMGs 
and their various subclasses at all wavelengths and redshifts 
and in various environments. Even their star formation rates 
(SFRs) remain uncertain in most cases because of a lack of 
direct observations at the FIR/submm wavelengths of their 
maximum emission. The identification of large samples of 
SMGs is important for investigating the properties of these 
galaxies (SFR, luminosity, spectral energy distribution [SED], 
stellar mass, AGN content, spatial structure, radio and X-ray 
parameters, clustering, etc.) on a statistical basis, as a function 
of their various subclasses, redshift, and environment. This is 
the main goal of the wide-field submm surveys planned with 
SCUBA2 and Herschel. 

Although Spitzer generally lacks the sensitivity to detect 
SMGs in the far-IR, its very good sensitivity in the mid-IR al- 
lows the efficient detection of a significant fraction of SMGs in 
the very large area observed by its wide-field surveys, and in par- 
ticular the ~ 49 deg^ Spitzer Wide-are a Infrared Extragalactic 
(SWIRE) survey (ILonsdale et al. 20031) . From an analysis of a 
sample of ~ 100 SMGs observed with Spitzer, Lo09 have es- 
timated that SWIRE has detected more than 180 SMGs with 
fi.2mm > 2.5 mJy per square degree at lApm and in several 
IRAC bands from 3.6 to 8.0jum. However, the identification of 
SMGs among SWIRE sources is not straightforward, since it re- 
quires inferring FIR emission from mid-IR photometry in ob- 
jects with various SEDs, especially as regards AGN versus star- 
burst emission, and various redshifts. 

We have therefore undertaken a systematic study of the 
1.2 mm emission from the best SMG candidates among Spitzer 
bright 24 yum sources, focusing on z ~ 2 s tarburst candidates. In 
Lonsdale et al. (20()6l) and Lo09 (see also Weedman et al. 2005 
Farrah et al. 20081) . It is shown that selecting sources with a sec- 
ondary maximum emission in one of the intermediate IRAC 
bands at 4.5 or 5.8yum provides an efficient discrimination 
against AGN power-law SEDs. In particular, 24 /urn bright 
"5.8 yum-peakers" have a high probability of being dominated by 
a strong starburst at z ~ 2, whose intense 1 .1 pm feature is red- 
shifted into the 24 pm band. A first 1 .2 mm MAMBO study of a 
sample of ~ 60 bright SWIRE sources (Lo09) has confirmed that 
such a selection yields a high detection rate at 1 .2 mm and a sig- 
nificant average 1.2 mm flux density, showing that the majority 



of such sources are high-z ULIRGs, probably at z ~ 2. However, 
as described in Lo09, this sample was selected with the aim of 
trying to observe the "5.8yum-peakers" with the strongest mm 
flux over more than 10 deg^. This was achieved by deriving pho- 
tometric redshifts, estimating the expected 1.2 mm flux densities 
by fitting templates of various local starbursts and ULIRGs to the 
optical and infrared (3.6-24yum) bands, and selecting the candi- 
dates predicted to give the strongest mm emission. Therefore, the 
selection criteria of this sample were biased, especially toward 
the strongest 24pm sources and those in clean environments. 
We report here the results of an analogous MAMBO study, but 
of a complete 24yum-flux limited sample of all SWIRE "5.8 pm- 
peakers" in a 0.5 deg^ region within the SWIRE Lockman Hole 
field, with F24^m > 400juJy and r'y^^^ > 23 (see Sec. 2 for 
a precise definition of "5.8jum-peakers", which of course de- 
pends on the actual SWIRE data and limits of sensitivity and 
accuracy). This region was selected because of the richness in 
multi-wavelength dat a, in particular the exceptionally deep radio 
data at 20 cm (VLA, lOwen & Mortso n 2008) , 50 cm (G MRT, 
Owen et al. in prep.), and 90 cm (VLA, lOwen et al. 20091) . Our 
study aims at characterizing the average multi-wavelength prop- 
erties of these sources, their dominant emission processes (star- 
burst or AGN), their stellar masses, and their star formation 
rates. We adopt a standar d flat cosmology: H q-TI km s ' Mpc 
Qm=0.27 and Qa=0.73 dSpergel et al. 2003h . 

2. Sample selection and ancillary data 

We selected all Spitzer/SWIRE "5.8jum-peakers" with 
/^24^m > 400//Jy in the 42' x 42' (0.49 deg^) J1046+59 field 
in the SWIRE Lockman Hole, centered at Q'2000 - 10h46m00s, 
^2000 - +59°0r00" (Fig[T]i- A source is considered to be a 
"5.8 jum-peaker" if it satisfies the following conditions: e^im < 
F4.5iim < -Fg.g^tn > without Consideration of uncer- 

tainties. 13 sources have no detection in the 8.0//m band. We 
assume that these sources are also "5.8//m-peakers" because 
their fluxes at 5.8jum are greater than the detection limits at 
8.0jum (<40juJy). We also require that the sources are opti- 
cally faint, i.e., > 23, to remove low redshift interlopers 
d Lonsdaleetal.2()() ^) . The selected sample contains 33 sources, 
which represent 6% of all sources with ^24^^11, > 400 //Jy and 
r' > 23 in the field. 

Vega 

We use the SWIRE internal catalogue available at the time of 
definition of the project, September 20 06. Details on the SW IRE 
observations and data are available in lSurace et al. (2()05ft . The 
2006 SWIRE internal catalogue has been superseded by the cur- 
rent version, dated 2007. Since our selection, MAMBO observa- 
tions, and analysis were based on the 2006 catalogue, but there 
are no significant changes for our sources in the 2007 catalogue, 
and the numbers of sources selected in the two versions of the 
catalogue do not vary significantly, for this work we will use the 
selection from 2006 data. Based on the analysis of the sources 
that would have been missed or included applying our selection 
criteria to the two versions of the catalogues, we find that the 
original sample selected from the 2006 SWIRE catalogue re- 
mains representative of a sample strictly meeting our selection 
criteria. 

We have reported in Table[T]the fluxes from the 2006 SWIRE 
catalogue. The optical magnitudes have been obtained with the 
MOSAIC camera on the 4-meter Mayall Telescope at Kitt Peak 
National Observatory (e. g., Mulleretal. 1998). However, for 
a few sources, the optical data diff'er slightly from those avail- 
able in the SWIRE catalogue because a measurement at each 
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Fig. 1. Positions of the 33 sources of our sample (42' x 
42' J 1046+59 field, centred at the VLA position observed 
by lOwen & Morrison ( 2008). Q'2000 — 10h46m00s, ^2000 — 
+59°OrOO")- The red filled triangles are the 4cr detections at 
1.2 mm. The green filled circles are the 3cr tentative detec- 
tions at 1.2 mm. The blue squares are sources with a signal at 
1.2 mm lower than 3cr. The sizes of symbols are proportional to 
1.2 mm strengths. The op en symbols show no detection at 20 cm 
dOwen & Morrison 20081) . 



source position was performed for all non-detections in the cata- 
logue. The revised optical magnitudes are listed in Table [1] The 
imaging in the central 0.3 deg^ is much deep er (by ~ 0.7 mag). 
Exact depths are listed in Pollett a et al. (20061) . In addition to the 
SWIRE data, we obtained near-i nfrared (NIR) data from two 
sources: 1) from Palomar/WIRC jWilson et al. 2003h : 2) from 
UKIRT/WFCAM jHenrv et al. 20001) in the J (A^ 1.25 jum), H 
{A - 1.63 jum), and K {A- 2.20jt/m) bands as part of a NIR sur- 
vey of the field. A description of these data and their reduction 
will be published in Strazzullo et al., in prep. The NIR data are 
reported in Table [1] 23 sources are detected in all three bands. 
24 sources are detected at least in the jiT- band. The WFCAM K 
band data are public data from UKIDSS ( Lawrence et al. 2007t 
[Warren et al. 2007) . Pipeline process ing and the science archive 
are described bv lHamblv et al. (2008i) and Irwin et al. (in prep.). 

Table |2] compares the selection criteria for our sample 
to those for similar samples of bright 24jum sources. The 
sample of Lo09 is based on the same criteria, but is bi- 
ased toward sources brighter at 24 /urn, with 819/zJy on av- 
erage vs. 566 ;uJy for the present sample. The sample of 
iFarrah et al. (2008*) is similar but aim ed at "4.5 //m-peak ers"; the 



sample of Huang et al. (2009) and Yo unger et al. (20091) is based 
on different IRAC criteria, but they indeed select almost exclu- 
sively "5.8 yum -peakers" (Sec. 5.1). On th e oth er hand, the selec - 
tion criteria of lMagliocchetti et al. (20071) and lYan et al. (20051) . 
which do not use the IRAC flux densities, do not discriminate 
against AGN and yield a large proportion of AGN. 

Compared to the twin starburst sample of Lo09, the present 
sample is complete down to a 24 yum flux density of 400 yuJy. It 



thus includes weaker 24 jt/m sources on average, but it should be 
free from the selection biases present in the Lo09 sample that 
resulted from the effort to maximize the number of detections at 
1 .2 mm. In addition, our sample benefits from very deep radio 
data at 1400, 610, and 324 MHz (Sec 4.3). The positions of the 
sources in the field are shown in Figure [T] The radio flux densi- 
ties are listed in Tabled 



3. MAMBO observations and results 

Observations were made during the winter 2006/2007 
MAMBO observing pool between December 2006 and March 
2007 at the IRAM 30 m telescope, located at Pico Veleta, 
Spain, using the 117-element version of the MAMBO aiTay 
dKrevsa et al. 1998) operating at 1.2 mm (250 GHz). We used a 
standard "on-off" photometry observing mode with a secondary 
mirror wobbling at a frequency of 0.5 Hz between the source and 
a blank sky position offset in azimuth by ~ 35". Periodically, 
the telescope was nodded so that the sky position lay on the other 
side of the source position. Pointing and focus were updated reg- 
ularly on standard sources. Nearly every hour, the atmospheric 
opacity was measured by observing the sky at six elevations. 
Our observations are divided in 16 or 20 "subscans" of 60 sec- 
onds each. In this operating mode, the integration time is ~ 40 s 
(20 s on source and 20 s on sky) per subscan. Observations of 
each source were never concentrated in a single night, but dis- 
tributed over several nights in order to reduce the risks of system- 
atic effects. The initial aim was to observe the 32 sources (L-12 
was observed in the project described by Lo09, under the name 
"LH-H") with an rms ~0.8mJy, which corresponds to ~0.6h 
of integration for the system sensitivity ~ 35-40 mJy s"^^^ in av- 
erage weather conditions. As seen in Table H] this was achieved 
for most of the sources. However, for ~ 20% of the sources the 
rms was instead ~ 0.9- 1 . 1 mJy, while a similar number of sources 
were observed somewhat longer to reach an rms ~ 0.5-0.6 mJy 
in order to confirm their detection. 

The data reduction is straightforward thanks to the MOPSIC 
software packag^l which is regularly updated on the MAMBO 
pool page. This program reduces the noise due to the sky emis- 
sion if it is sufficiently correlated between the different bolome- 
ters. This process is generally sufficient for the majority of ob- 
servations. However, in a few cases, some scans may present 
faults due, e.g., to lack of helium in the cryostat or problems 
of acquisition. These scans are validated or rejected after close 
inspection. 

The results of our observations at 1.2 mm are reported in 
Table |4l The average flux density (with equal weight) of the 
entire sample is 1.5 6 + 0.22 mJy, very c omparable to Lo09 
(1.49 + 0.18 mJy) and 'Younger et al. (2009) (1.6 + 0.1 mJy), but 
greater than obtained by Lutz et al. (2005) for a Spitzer-selected 
sample of high-z starbursts and (mostly) AGNs (0.74+0.09 mJy). 
The median for our sample is 1.61 mJy. This confirms that on 
average the majority of these sources are SMGs (at z ~ 2 , 
1.6 mJy coiTesponds to ~4mJy at 850//m: iGreve et al. 20041) . 
Because of the limited integration time, only four sources were 
solidly detected at >4cr. However, the fraction of sources at 
least tentatively detected at > 3cr is 39% (13/33 sources). It is 
worth stressing that the reliability of such 3cr tentative detec- 
tions in careful "on-off" MAMBO observations is much higher 
than those detected in a mm/submm map amongst hundreds of 
possible resolution elements, where flux boosting is inevitable. 



Documentation 
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This on-off 3cr 'detection' rate is similar to what was ob- 
tai ned by Lo09 for a similar sample, higher than obtained 
by iLutz et al. (2005h for their 5/?/fzer-selected sample of star- 
bu rsts and (mostl y) AGNs (18%), and lower than obtained 
by lYounger et al. (2009) with deeper observation of a similar 
sample of Spitzer-selected z ~2 stai-bursts (75%). 

4. Source properties 

4.1. Spectral energy distributions and redshifts 



is based on spectroscopic redshifts (Fig.|2]i. Thus, all these works 
select sources in a similar redshift range. Indeed, the actual 
redshift distributio n of our sample might be similar to that of 
iHuang et al. (2009h and lYounger et al. (20091) and concentrated 
within a rather narrow redshift range around z ~ 1.7-2.3 as 
shown by the few spectroscopic redshifts reported for the sample 
ofLo09. 

Twelve sources from our sample have redshifts 
from the catalogue of SW IRE photometric redshifts of 
iRowan-Robinson et al. (20081 ). For ten of these sources, our 
photometric redshifts determined by the SED fitting show 
a good (+ 10%) agreement with the determinations from 
iRowan-Robinson et al. (20081) . 
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Fig. 2. Histogram of redshifts for our sample (photometric, 
thick solid red line), the full sample from Lo09 (photomet- 
ric ar id spectroscopic, do tted black line), and the sample 
from lYounger et al. (20091) (spectroscopic, long-short-dashed 
blue line). 



In order to estimate photometric redshifts, we fit the spectral 
energy distributions (SEDs), including optical and IR (< 24yum) 
data for each source, with a library of galaxy templates follow- 
ing the method described in Lo09 and Polletta et al. (20 0j). Th e 
SEDs are fitted using the Hyper-z code (B olzonella et al. 2000h . 
and the effects of dust extinction are taken into account. As dis- 
cussed in Lo09, such photometric redshifts are limited in accu- 
racy and have uncertainties of +0.5. 

SED fits for all the sources in the sample are shown in the 
Appendix (only in the electronic edition), and the photometric 
redshifts are listed in Table |4] The photometric redshift distri- 
bution of the sample is shown in Fig. |2] All redshifts but four 
are within the range 1.5 < z < 2.5. The average redshift from 
these SEDs is < z >= 2.08 (median = 2.04, scatter = 0.32, and 
semi-inter-quartile range = 0.26). This result is consistent with 
our selection criteria, which assume that the 7.7 yum PAH band is 
redshifted into the 24 /im MIPS band and the 1 .6 fim stellar band 
into the 5.8 fim IRAC band. The redshift distribution of our sam- 
ple is similar to the one measured in Lo09 (< z >= 1 .97 + 0.05), 
which is mo stly based on photo metric redshifts, and the one 
measured in lYounger et al. (20091) (< z >- 1.96 + 0.10), which 



4.2. Comparison between 1.2mm and 24 fim flux densities 
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Fig. 3. Observed MAMBO 1 .2 mm flux density as a function of 
24 yum flux density. The black squares show the average flux den- 
sities for the first F24//m quartile (8 sources) and for the 25 other 
sources, respectively. The filled (open) red circles are the sources 
that are 4cr detections (3cr tentative detections) at 1.2 mm. The 
green arrows are 2 cr upper limits for the sources that are not 
detected at 1.2 mm. 



In order to investigate whether there is a coiTelation be- 
tween mid-IR and mm emission, we compare the flux densities 
at 1 .2 mm and at 24 fim (Fi.2mm and F24f2m) in Figure[3] Because 
of the limited sensitivity of the 1 .2 mm data, we have stacked the 
data for the first F24/jm quartile (8 sources), and independently 
for the 25 other sources. These stacked values have been com- 
puted with the observed /^i.2mm values. Fig. |3] shows that it is 
impossible to see whether there is a coiTelation between mid-IR 
and mm emission with this sensitivity. However, the average for 
the highest F24/jm quartile seems ~ 1 .5 - 2 times larger than the 
average for all the other sources. 

Because of the 24yum selection, the /^i.2mm/^'24;/m ratio of 
our sources is relatively low compared to that submm selected 
SMGs, as in the case of Lo09. The ratio of the average Fi.2mm 
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to the average F24/im is 2.76 + 0.50 for the entire sample of 33 
sources (and 4.94 + 0.51 for 13 sources with 1.2nun S/N > 3, see 
Table 3). As seen in Fig. [8] these ratios of averages are a fac- 
tor ~ 4 (~ 2) smaller than that of a sample of literature SMGs 
(Sec. 5.2), and a factor ~ 6 (~ 10) larger thai i one of the (AGN 
dominated) sample of bright 24 /vm sources of lLutz et al. (20051) . 
They are a factor ~ 10 (~ 5) smaller than that of the extreme 
SMG template Arp 220 at z = 2, but comparable to those of star- 
burst templates (M 82 or NGC 6090) and AGN-starburst com- 
posites (IRAS 19254-7245 South). 

4.3. Radio properties and tlie nature of tiie sources 




100 200 300 

^20om (Wy) 

Fig. 4. Observed MAMBO 1.2 mm flux density as a function 
of 20 cm flux density. The large black symbols show different 
stacked values: the entire sample (filled triangle), all sources 
with > 3cr (filled diamond), the sources with > 3cr signal 
without the two strongest radio sources (open square), and the 
sources with < 3 o" signal (filled square). The filled (open) red 
circles are the sources that are 4cr detections (3cr tentative detec- 
tions) at 1 .2 mm. The green arrows are 2 cr upper limits for the 
sources with S/N < 3 at 1 .2 mm. The sohd black line shows the 
correlation derived from IChapman et al. (20051) for z - 2 (as- 
suming Ffi,50ij„i/Fi,2mm - 2.5 atz~ 2). The dashed lines are the 
limits of this correlation at +lcr. The labelled sources are those 
with the highest 20 cm/1.2 mm flux density ratios (see Table|4]i. 

The studied sources benefit from exceptionally 
deep VLA data at 1.4 GHz (rms = 2.7//Jy in the cen- 
ter of th e field, 12-15/L(Jy in most of the 0.5 deg^ field; 
lOwen & Morrison 2008). Such a depth yields radio detections 
for almost the entire sample. 8 of 33 sources are not detected due 
to a loss of radio sensitivity in the outer parts of the field, largely 
due to the decrease in primary beam sensitivity and bandwidth 
smearing. The GMRT 6 10 MHz observations are also very 
deep (rms = lOyuJy), cover almost the entire field, and detect 
all 33 sources. The VLA 324 MHz observa tions reach a dept h 
of rms = 70/iJy and cover the entire field jOwen et al. 20091) . 



They yield detections for 17 sources. The radio flux densities 
are listed in Table |4] The cross identification was made by 
comparing the radio and the SW IRE positions. Followin g 
the method of llvison et al. (20071) and (Pownes et al. (19861) . 
we have verified that all our sources have very reliable radio 
associations, with an averag e probability of spu rious association 
in 2" of <f> = 0.001. For llvison et al. (2007h . the association 
isreUable if P < 0.05. 

Based on the coiTelation between Lfir and the radio 
luminosi ty in star-forming regions and in local starburst 



galaxies jHelou et al. 1985HCondon 1992l;ICrawford et al. 19961; 
ISanders & Mirabel 1996 ^). the radio and FIR luminositie s are 
expected to be linked as well at z ~ 2 jlbar et al. 20081) . We 
discuss this in Sec. 4.4, together with the derivation of Lfir- 
However, if this correlation is verified at z ~ 2, one may ex- 
pect a straightforward relation between the radio and 1.2 mm 
flux densities for starbursts with similar SEDs. To check that, 
we plot the 1 .2 mm flux density as a function of the radio flux 
density at 20 cm in Fig. |4] Because of the limited sensitivity 
at 1.2 mm, we consider average values of 1.2 mm flux density 
and radio flux density, for the entire sample, for the 13 sources 
with 1 .2 mm S/N > 3, and for the 20 sources with 1 .2 mm S/N < 3 
(Table [3]l. In Fig. |4] we also report the ratio of the a verage flux 
densities and Icr dispersion found by [C hapman et al. (2005h (see 
alsoe^ Condon 19921; [Smail et al. 20021; lYun & CarilU 20021; 
llvison et al. 2002) for a sample of radio-detected sub-millimeter 
galaxies at z ~ 2. The mm/radio ratios for both the whole sample 
and the S/N > 3 sources are reasonably compatible with the cor- 
relation bet ween mm/submm an d radio fluxes found for z ~ 2 
starbursts (IChapman et al. 2005h . However, the ratio between 
average values of 1 .2 mm flux density and radio flux density, for 
the entire sample and especially for the 20 sources with 1 .2mm 
S/N < 3, might be slightly smaller compared to submm selected 
galaxies. This could be explained as an effect of the 24yum selec- 
tion, and by a greater AGN contribution and/or hotter dust (see 
Sec. 5.3 and 5.4 for a complete discussion). 

In order to assess whether the radio emission observed in 
our sources is associated with AGN or with star-forming re- 
gions, one may also consider the radio spectral shape, the rest- 
frame radio luminos i ty, and the radio m orphology (see e.g. 
iBiggs & Ivison 20"08t ISevmour et al. 2008b . The radio spectral 
index a, where Fy oc V, is first calculated from the flux densi- 
ties at 20 cm, 50 cm, and 90 cm, when they are available, using 
the best power law fit between these three wavelengths, and is re- 
ported in Table |4| When only two radio fluxes are available, the 
index is just derived from the flux density ratio. Seven sources 
have no determined spectral radio index because of a lack of ra- 
dio detections at 1.4 GHz and 324 MHz. 

Most of our sources have a radio spectral index in the 
range ~ -0.4 to -1.2 (average < a >- -0.64 + 0.07; 
median - -0.74). This is not very discriminating since, con- 
sidering the uncertainties, such values of a are typical for 
star-form ing galaxies, t ype II AGN, and rnany radio galax- 
iesJe^'Condon 19921 'Po Uetta et al. 2000l: ICihegi et al. 20031: 
iDe Breuc k et al. 2000, 2001). About 20% of flie sources, mostly 
among those with some radio excess, could either be in the same 
range or have a > -0.5 typical of flat spectrum sources. 

The radio luminosity and the FIR-radio correlation are most 
often expressed in term of the rest-frame luminosity at 1 .4 GHz, 
L 1.4 GHz- With the sa me assumptions as for the SMG sample of 
iKovacs et al. (20061) . we can rewrite their Eq. (7) as 



il.4GHz(WHz-l)=47rD2F20c 



(-a-\) 



(1) 
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where Di is the luminosity distance. In cases with detections at 
at least two radio wavelengths, we assume the derived value of a. 
When only Fsocm is available, we assume a fixed a value equal 
to the average value found for our sample, < a >= -0.64 
Using the flux at 6 10 MHz, this equation becomes 



Li.4GHz(WHz-i) =4;rDiF50c 



1 



2.3 



I +z I +z 



(2) 



Finally, we examine the radio sizes of our sources. A 
large size (> 1") may indicate radio jets or lobes asso- 
ciated with the presence of a radio loud AGN. However, 
it is also possible that the extended radio emission is as- 
sociated with star-forming regions distributed across the 
galaxy, as is the case for some local ULIRGs (e.g., 
iMurphv et al. 200 ij) and a large fraction of classical SMGs, e.g. . 



-70% in 'Chapman et al. (2004) (see also [Seymour et al. 20081 
[Richards et al. 2007: Biggs & Ivison 2008 ). Nine sources are 



extended (major axis >5kpc) in the VLA 20 cm images, with 
most sizes > lOkpc, and the remainder are either unresolved or 
not detected (see Table|4]l. The radio luminosities of the nine ex- 
tended sources range from ~ 10^^'^ W Hz"' to ~ lO^"*'^ W Hz"' . 



4.4. Far-infrared luminosity 

Estimates of Lfir for our individual sources are quite uncer- 
tain due to the lack of data between 24 jjm and 1 .2 mm, where 
most of the far-infrared energy is emitted. None of our sources 
are detected at 70;um or 160yL(m at the SWIRE sensitivities (the 
3cr limits at 70 and 160 yum are 18 and 108 mJy, respectively). 
There are not even detections in the ~ 2.5 times deeper MIPS 
images available from GO Program 30391 (PI: F. Owen) in the 
center of the field, which is well observed at 20 cm. We thus 
make use of these deep MIPS data only for stacking the 70 and 
160yum images to constrain the average Lfir of our sample. 

We estimate average flux densities from the stacked MIPS 
images at 24, 70, and 160yum, following the method described 
in Lo09. An initial stack based on the SWIRE MIPS images was 
made for the entire sample of 33 sources, but it yields only a 
marginally significant (~2-3 cr) detection at 160 jum. The aver- 
age flux densities at 70/im and 160jum are almost 5 times lower 
than the SWIRE 3 cr limits (Table |6]l. We thus use stacks made 
with the GO-30391 images of the 21 sources covered by these 
deeper MIPS observations. The results from stacks are reported 
in Table |6] In addition to the stacked flux densities for the entire 
sample, we also report stacked flux densities for two subsamples 
apiece: 10 sources with 1.2 mm SNR>3, and 11 sources with 
SNR<3. 

The median stacked flux densities at IQfim and 160yum for 
the entire sample are quite faint, about 3 mJy and 13 mJy respec- 
tively. We use these median stacked flux densities in the follow- 
ing analysis. There is no appreciable difference in < Fieo/jm > 
and < F70;jm> between sources tentatively detected (>3 cr) and 
undetected (<3cr) at 1.2 mm. The median 160 jum and 70yum 
flu x densities seem sma ller than for the previous similar sample 
of lYounger e t al. (20091 ), which has higher average 24 //m flux 
density (Table |6ll. 

We combine the median flux densities at lOjjm, 160;um, 
and 1.2 mm to build a median far-IR SED for our sample. We 
fit the median FIR SED assuming the average redshift of the 
sample (z = 2.08) with a "graybody" model with a fixed value 
of the emissivity i ndex ( f5 - 1.5, see e.g. Kovacs et al. 2006; 
iBeelen et al. 20061) to derive the temperature, Tjust, of a sin- 
gle dust component. The best-fit value for the whole sample 



is Ldust = 37 ± 8K. This model yields Lfir = 2.5 x IO'^Lq, 
which is well in the ULIRG range and may be considered the 
average value for "5.8yL(m-peakers" with F24/jm > 400jt/Jy and 
r > 23. Note that these values for Tdust and < Lfir > are slightly 
smaller than the valu es of Lo09, Tdust = 4 1 K and < Lfir >= 
4.6 X lO'^Lo and of lYounger et al. (20091) . Tdust = 41 K and 
< Lfir >= 3.8 X 1012 Lq. 

The actual Lfir of our sources will span some range around 
this value. Although it is not possible to derive accurate values 
of Lfir for each source in the absence of individual detections at 
160jum or at another wavelength close to the FIR maximum of 
the SED, several alternative approaches are possible for estimat- 
ing individual Lfir. 

One can simply infer Lfir from the flux density at 1.2 mm, 
Fi.2mm- The value of Fi.2mm can place strong constraints on 
Lfir if we assume that the whole Lfir and the radiation de- 
tected at 1 .2 mm are produced by dust heated by the same mech- 
anism. The relation between Lfir and 2 mm can be derived as- 
suming a thermal spectral model. A "graybody" model, with 
a single dust component and emissivity index, e.g. fi - 1.5 
(Beelen et al. 2006), yields: 



LFlRiLe) ^KX 10'2Fi.2mm(mJy) 



(3) 



where the proportionality factor, a- ~ 1 - 5, depends mainly on 
Tdust but also by a factor of a few on redshift. For redshifts ~ 2, 
K ^ 1.3 for Ldust = 35K and = 2.14 for Ldust = 40 K. The 
dust temperature, Ldust, is difficult to constrain without obser- 
vations in the FIR-submm range, e.g., at 160 or 350;um. Ldust 
is known to spa n a large range, fro m 21 K to 60 K in 350 /im- 
detected SMGs (iKovacs et al. 20061) . and from 34 K to 47 K in 
z ~ 2 Spitzer-select ed and mm-detected star-forming galaxies 
(lYounger et al. 20091) . In the absence of other information, we 
assume the average values inferred from our stacks, Ldust = 
37 + 8 K for the whole sample, and 36 + 3 K and 39 + 2 K for the 
samples with 1 .2 mm S/N > 3 and < 3, respectively. Assuming a 
"graybody" model with Ldust - 37 K, 36 K and 39 K, the values 
of the factor k are 1.57, 1.41, and 1.93, respectively. Lfir values 
thus derived using Eq.|3] LFiR mm, are reported in Table |5]for the 
1 .2 mm 4cr detections and 3 cr tentative detections. For the other 
non-detections, we report upper limits to Lfir. 

We also estimate Lfir of each source using the avail- 
able radio data and the well-k nown FIR-radio relation for 
local starbursts ( Condon 19921) . The original definition is 
based on the 60ju m and 100/zm flux densities, but following 
ISaiina et al. (20081) . we have adopted the definition: 



^ = log 



120 



log 14.03 
^ IWHz-i/ 



(4) 



where L4o_i2o is the integrated IR luminosity between 40 pm and 
120 jum in the rest frame. We assume LFiR,mm = i'40-i20- 

Although Eq. (4) has been found to hold for local sources, 
its validity se ems to be confirmed al so at the re dshifts of our 
sources (Kov acs et al. 2 006: Younger et al. 20091) . The relation 
between Lfir and the 20 cm or 50 cm flux density is thus straight- 
forward, if we assume the typical value q = 2.34 and LFiR,mm = 

^^40-120 



Lfir(Lq) 



2.04X 10-'2Li.4ghz(WHz-1). 



(5) 



However, this method is valid only if the radio emission is 
not significantly affected by an AGN. This expression will give 
an upper limit to the star-formation related Lfir if an AGN con- 
tributes significantly to the radio emission. The Lfir derived us- 
ing Eq.|5] LFiR,radio, are reported in Table |5] 
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Fig. 5. The radio-FIR q factor (FIR-to-radio ratio) (Eq. 4). The 
large black symbols show the stacked values of different sam- 
ples: the entire sample (33 sources, filled circle), the 13 sources 
with S/N > 3 at 1.2 mm (filled square), and the 20 sources with 
S/N < 3 (open square). The solid black line shows the typical 
value q - 2.3 4 for local starbursts. The dotted b lack lines are 
the 3cr spread (lYun et al. 20011; ISaiina et al. 20081) . q values for 
individual sources are also shown for the sources with > Act de- 
tections at 1 .2 mm (small filled symbols) and the sources with 
> 3cr signal (small open symbols). The red triangles show the 
sources identified as having the highest 20 cm/1. 2 mm flux den- 
sity ratios. The red circles are the sources with no 20 cm excess. 



The values of these two estimates of LpiR are within a factor 
of 3 for most of the sample. They range from ~ 0.5 to ~ 10 x 
10'^ Lq, confirming that almost all our sources are ULIRGs with 
luminosities greater than 10'^ Lq. The average Lfir derived from 
the radio flux densities is < LFiR.radio >= (4.14 + 0.53) x 10'^ Lq. 
Since this value may be slightly overestimated because of AGN, 
it is consistent with LpiR derived from fitting the median FIR 
SED, < Lfir >= 2.5 x lO'^L^. 

Alternatively, the radio-FIR relation (Eq. 5) can be applied to 
identify or test for the presence of AGN-driven radio activity in 
our sample (similarly to Sec. 4.3 but slightly more rigorously). 
The agreement or deviation from this correlation can be easily 
expressed by the value of the ^-factor, as reported in Table |5] 
We obtain an average value < q >- 2.11+0.12 for the whole 
sample. 

This is not very different from the av erage value < q >= 
2.34 found for local star-forming galaxies jYun et al. 200 lb . As 
shown in Fig |5] the <7-factors of the three sources identified as 
having the highest 20 cm/1. 2 mm flux density ratios seem to be, 
on average, lower than those of most other sources. This is not 
surprising, since the q factor is another way to quantify the radio 
excess. 



4.5. Stellar mass 

The IRAC-based selection of our sample corresponds to 
a rest-frame NIR selection; therefore, if we ignore a possi- 
bly significant AGN contribution at these wavelengths, we are 
directly sampling the stellar component (Lo09). Applying the 
method developed by Berta et al. (20041) . we have derived an 
estimation of the stellar mass for our sources. As summa- 
rized in Lo09, this derivation uses not only the IRAC data, 
but also optical and 24jum flux densities to perform a mixed 
stellar population spectro-photometric synthesis, and to con- 
strain t he extinction and the luminosity of young stellar popu- 
lations jBerta et al. 20041) . A redshift - here photometric - is 
also required. The stellar m asses, initially computed with a 
0.15 - 120MM lSalpeter (1955) initial mass function (IMF), are 
converted to a .10 - 100 Mq I MF by applying a correction 
factor of 1.162 (iBerta et al. 2003i) . All the stellar masses are 
shown in Table|5] Their values range from 0.8 to 7.2x10" Mq, 
with a median value of 1.37x10" Mq, and an average value 
1.77x10" Mq. 

As discussed in Lo09, these values may be somewhat 
overestimated for several reasons: our models may not fully 
take into accou nt the TP-AGB contribution to infrared light 
(Marasto n 20051) ; they assume a Salpeter-like IMF, which yields 
higher masses than other IMFs; and although the near-IR SED is 
dominated by stellar emission, it is possible that an AGN com- 
ponent is also present. Therefore, these estimates should be con- 
sidered as upper limits to the true stellar masses, and may be 
overestimated by a small factor up to ~2. 

Using the same method, Lo09 found a median stellar mass 
- 1.80x10" and an average value 2.16x10" Mq, about 
30% larger than for our sample. This is consistent with the 
comparison of the average values of the 5.8 /im flux density of 
both samples, which reflects the rest-frame luminosity at 1 .6 fim: 
53 yuJy for our sample and 77 //Jy for that of Lo09. Indeed, as 
discussed by Lo09, the direct comparison of the NIR rest-frame 
luminosity at 1.6/^m is probably the most consistent way to 
compare stellar masses in different samples, in order to avoid 
the dependence on the various methods applied to infer stellar 
masses from infrared fluxes. We have thus estimated the rest- 
frame luminosities at 1.6;um, vLv(l-6yum) (Table|5]), by interpo- 
lating the observed IRAC fluxes as in Lo09. As expected, the 
average rest-frame luminosity at 1.6yum in Lo09, 4.1x10" Lq, 
is a factor of 1.6 larger than for our sample, 2.5x10" Lq (Table 
|5j. The sources in our sample are as luminous as, or slightly 
more luminous at 1.6yum by a factor ~ 1.5 than, submm selected 
SMGs, and should also be ~ 1 .5 times more massive than clas- 
sical SMGs assuming the same mass-to-light ratio of Lo09. Our 
sample also shows a mass-to-light rati o consistent with the radio 
galaxy sample detected at 24 /urn by S eymour et al. (20071 ). who 
computed stellar masses directly from the luminosity at 1 .6 i-im 
in the rest-frame. 

4.6. Star formation rate 

From the estimated Lfir, we derive the star forma- 
tion rate (SFR) of our sources assuming the relation 
from lKennicutt(1998h using a 0.10 - 100 Mq Salpeter IMF: 

SFR(Mo yr-^) « 1.8 x 10""'Lfir(Lo). (6) 

The SFR values derived from LFiR,mm and LpiR^iadio are listed 
in Table |5] The SFRs derived from LFiR,iadio for the sources with 
possible radio excess may be overestimated due to the possi- 
ble AGN contribution to their radio emission. The mean SFR 
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Fig. 6. {Adapted from Fig. 14 o APaddi et al. 2007\} . Star forma- 
tion rate vs. stellar mass. Values of SFR deduced from radio data 
for the sources of our sample are represented by the red trian- 
gles (see Sees. 4.4, 4.6 and Table|5]l. The aiTows are 2cr limits. 
The large black symbols show stacked values for different sam- 
ples: all sources (radio determination = filled triangle, millimeter 
determination = open triangle), the sources with >3 cr signal at 
1.2 mm (radio determination = filled diamond, millimeter deter- 
mination - open diamond), and the sources with <3 cr signal at 
1 .2 mm (radio determination - filled circle, millimeter determi- 
nation = open circle). The small bla ck squares ar e for the 24yum 
BzK sources from GOODS (Daddi et al. 20071) . The magenta 
star shows a typical value for SMGs. The large green squares 
trace the average SFR-mass relation in GOODS-N (160 arcmin^) 
determined from radio stacking of K<2Q.5 galaxies in three 
mass bins; the blue Une is SFR = 200 M^-^ (Mq yr~'), wh ere Mu 
is the stellar mass in units of 10" Mq jDaddi et al. 2007h . 



derived by averaging the radio-based estimates of all sources 
is ~75OM0yr"'; the mean derived from the average LpiR is 
~450Moyr-i. 

To estimate the contribution of "5.8/im-peakers" with 
>400yuJy to the star formation rate density (SFRD) 
of the universe, we consider the space density and aver- 
age SFR of our sources. Assuming a redshift interval of 
1.5< z <2.5, and an average SFR of 45OM0yr"', we de- 
rive the contribution of our sample to the (comoving) SFRD 
to be ~1.5-4xlO"^M0 yr"' Mpc"^. This value corresponds to 
-5% of the SFRD of all classical SMGs dAretxaga et al. 20071: 
IChapman et al. 20051) . i.e., to -10% of the SFRD of SMGs in 
the interval 1.5< z <2.5, and up to ~15% in the interval of 
1.7< z <2.3 where "5.8/zm-peakers" are mostly confined. 

The specific star formation rate (SSFR), defined as SFR/M*, 
ranges from ~10"'*yr"' to ~10"^yr"' for our sample. There 
is a tendency for the sources with the lowest stellar masses 
to have the highest SSFRs. This result is consistent with, e.g., 
[Noeske et al. (2007). However, it is seen in Fig.|6]that "5.8 fim- 
peakers" have significantly higher SSFRs than classical sBzK 
galaxies of comparable masses on average (Sec. 5.3). 



5.1. Comparison witii otiner sampies: 1) Spitzer selectioti 
(see Table\B 
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Fig. 7. IRA C color-color diag ram f or several samples. The 
sample of iHuang et al. (20091) and lYounger et al. (20091) is 
represented by the solid blue triangles. The red circles 
show our sample. T he open green stars are the sample of 
iLonsdale et al. (20091) . The arrows are 3cr upper limits for the 
di fference [3.6]-[8.0]. The inset box shows the selection criteria 
of lHuang et al. (2()09l) . Magnitudes are on the AB scale. 



The present work is an extension of the study presented 
in Lo09 on a sample of 61 sources selected in four SWIRE 
fields. The sources in Lo09 meet the same selection criteria 
and span the same redshift range, 1.4 < Zphot < 2.7, as our sam- 
ple. However, their selection is biased toward brighter sources 
at 24 /im, while the present sample is complete to a 24/im 
flux density limit of 400 /iJy. This bias results in an average 
24/im flux density a factor 1.5 greater than in our sample, 
< ^24//m > = 819yL(Jy vs 566yL(Jy (see § |2]i. Nevertheless, the 
1 .2 mm properties of the two samples are comparable. The aver- 
age flux densities at 1.2 mm are similar: 1.56+0. 22 mJy for our 
entire sample, against 1.49+0. 18 mJy for the sample in Lo09. 

Although the source selection in Lo09 was optimized to 
favor 1.2 mm bright sources, their detection rate is slightly 
lower than what we achieve with our complete sample, only 
26% (31% in their best observed field with a sensitivity sim- 
ilar to ours) with S/N>3 cr compared to 39% in our sam- 
ple. Both samples are characterized by a similar average 
LpiR determined using a thermal spectral model (see Sec 
4.4 and Eq.3), <LFiR,inm > = 2.5+0.4x10'^ Lq in this work vs 
(2.8±O.l)xlO'^L0 in Lo09 (restricted to the Lockman Hole 
field). Thus, in spite of the refinement implemented by Lo09 
in their sample selection to increase the chance of finding mm 
bright sources among "5.8//m peakers", the average mm pro- 
perties of their sample and our complete sample are consistent. 
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lYounger et al. (20091) have studied the FIR properties of a 
similar Spteer-selected sample, based on MIPS 70 and 160 //m 
detec tions and MAMB O 1.2 mm observations. This sample (see 
also iHuang et al. 2009b is selected based on IRAC colors, and 
24 /zm flux densities (see Fig.|2l Table |2] and §|2]i, yielding 12 
starburst galaxies with overall properties (redshifts, IRAC col- 
ors, 1.2 mm flux densities) similar to those in Lo09 and in this 
work. Consequently, the LpiR and derived star formation rates 
are also simil ar in all these sam ples. 

FinaUy, iLutz et al. (2005h studied the mm prop- 
erties of a Spitzer sample selected based on faint R 
band magnitudes, relatively bright 24/im flux densities 
(f 24u m > 1 mJy), and hig h F24^m/Ff,^m ratios jYan et al. 2005L 
120071: ISaiina et al. 20081) . This sample contains 40 sources 
at z ~2 exhibiting both starburst and AGN properties. 
They are on average fainter 1.2 mm emitters than our 
sources and they show significantly lower f'i.2mm/^^24/im 
flux density ratios, with < Fi 2mm >= 0.63+0. 20mJy (sim- 



100 



<F 



1 .2 mm 



>/<F24 



-- 0.43 (simple 



pie average) and 
averages), compared to < Fi. 2mm > = 1-56+0. 22mJy and 
< Fi,2mm >/< F24^im > = 2.76+0.50 for our sample (Table [3j see 
also Fig. [8]l. These differences are due to the fact that in the 
majority of these objects, the main source of power is an AGN 
rather than the powerful starburst required to produce bright 
mm flux, with the AGN dominating their 24 jum emission. 



5.2. Comparison witti otiier samples: 2) Submillimetre 
selection 

In order to characterize more quantitatively the difference in 
Fi.2mm/F24fim Aux density ratios between our sample of "5.8 fim- 
peakers" and classical SMGs, we show them as a function of 
redshift in Fig. |8] The sample of classical SMGs is the sub- 
set of the sample used in Lo09 completed by the SHADES 
SMG sample from lCoppin et al. (20061) and ilvison et al. (20071) 
detected at lAjim with 1.5<z<2.5. We also display a sub-set 
of "5.8 jum-peakers" from Lo09. Note that most of the redshifts 
available for the Lo09 and the SHADES samples are photomet- 
ric (Aretxagaet al. 2007). We also show in Fig. [8] the expected 
flux ratios for representative starburst and AGN templates. To 
facilitate the comparison, we also show the ratio of the aver- 
age flux densities for the samples of "5.8//m-peakers" from 
this work and from Lo09, and for a sample of classical SMGs 
(the SMG sample of Lo09 aug mented by those of SHADES: 
ICoppin et al. 20061 : llvison et al. 2 007). Our complete sample of 
"5.8/im-peakers" confirms that the ratio of the average 1.2 mm 
and 24 yum flux densities for this class of sources is smaller 
than for SMGs. However, the difference is slightly smaller than 
in Lo09 because of our lower mean 24yum flux density. The 
F\.2mmlF2Aiim Aux density ratios of our sample are also similar 
to one for the lensed LBG cB58, despite the fact that this lensed 
LBG is an or der of magnitude less luminous than the average of 
our sample ( Siana et al. 20081) . 

We have shown that "5.8//m-peakers" are z ~1 ULIRGs 
and that about 40% of them are bright mm sources. Thus, 
~ 40% of them also belong to the class of SMGs (Table 
|5]), one of the main classes of high-z ULIRGs. A detailed 
comparison between the MIR and FIR properties of a 
large sample of SM G s (iGreve et al. 2004 iPope et al. 2005 : 
Chapman et al. 20051: iBorvs et al. 20051: iFraver et al. 2004 



Hainline et al. 2009h and "5.8yum-peakers" is presented 



Lo09. In this study, Lo09 find that most "5.8 //m-peakers" repre- 
sent a sub-class of SMGs. The differences found by Lo09 with 
respect to classical SMGs are mainly related to the selection 
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Fig. 8. Observed flux density ratio, /^i.2mm/^24//m, versus red- 
shift. The small filled (open) magenta circles are the Act de- 
tections (tentative 3 cr detections) at 1.2 mm. The green arrows 
show the 2 cr upper- limits for the other non-detections. The 
crosses show values for a sample of SMGs (see text). 1.2 mm 
flux densities for this sample were derived from Fgso/jm when 
^^i.2mm is not available. The large symbols are the ratio of av- 
erage flux density vs average redshift for different samples: all 
33 of our sources (filled green square); all 13 of our sources 
with a > 3cr signal (magenta circle); the entire sample in Lo09 
restricted to the Lockman Hole field (open blue triangle); the 
sources with signal > 3crin Lo09 ( filled blue triangle ); the AGN- 
dominated 24 yum bright sample of iLutz et al. (20051) (black pen- 
tagon); and the SMG literature samp le (open star) (see Lo09; 
ICoppin et al. 2006t llvison et al. 20071) . Expected values for vari- 
ous starburst and AGN templates are also shown (Lo09): star- 
bursts Arp220 (solid green line), M82 (dot-dashed red line), 
and NGC 6090 (dotted magenta line); AGN-starburst composite 
IRAS 19254-7245 South (dashed black line); and AGN "QSOl" 
(short-long-dashed blue line). 



criteria for the "5.8//m-peakers". More specifically, the Spitzer 
selection favours sources with redshifts mostly concentrated in 
the range z ~ 1.7-2.3, rather higher stellar masses than classical 
SMGs at similar redshifts, brighter 24yum flux densities, and 
thus higher /^i.2mm/^^24/jm flux density ratios (see Fig. 9 in Lo09 
and Fig. [8j, and likely warmer dust temperatures. The mm- 
detected "5.8;um-peakers" and classical SMGs show similar 
mm/submm flux densities, implying that the main difference 
in the 1.2 mm/24 yum ratio comes from the 24 /im rather than 
the 1.2 mm intensity. This result implies that a high 24/im 
flux density does not directly translate into a high 1.2 mm flux 
density, and thus into high LpiR and SFR, in starburst galaxies at 
z ~2 (see also .Pope et al. 2008a; .Rieke et al. 2009.) . 



5.3. Nature of the sources: 1. Powerful starburst activity 

As discussed in Sec. 2, our selection criteria are devised to 
distinguish starburst ULIRGs from AGN through the presence 
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Fig. 9. Median optical through radio SED of all 33 "5.8 /vm- 
peakers" in our sample (red full circles). The 70yum and 160 yum 
median values correspond to the stacked values of the central 
21 sources (see § 14.4b . The magenta and cyan curves represent 
templates of IRAS 20551-4850 and Mrk231 normalized to the 
rest-frame luminosity at 1 .6 jum. The black dotted curve repre- 
sents a type 1 AGN template normalized at the lowest NIR dat a 
point (rest A = 2.4 //m) of the median SED dPoUetta et al. 20071) . 
The blue triangles represent the av erage fluxes of the sample of 
31 HDFN SMGs ( iPope et al. 2006h . The blue triangles at 70 //m 
and I60fj.m correspond to the stacke d fluxes fo r a subsample of 
26 sources from the same sample (Huvnh et al. 20 07). The green 
open squa res deUneate the av erage SED of the GOODS SMGs 
derived bv lPope et al. (2008bh . 



of the redshifted 1 .6 yum stellar "bump" which dominates over a 
strong AGN continu um in the near infrared. As shown by sim - 
ilar studies such as IWeedman et al. (20061) . iFarrah et al. (20081) . 
iHuang et al. (20091) . and Lo09, such criteria are very success- 
ful in selecting a majority of starbursts with strong mid-IR PAH 
features, and we have good evidence that this is also the case for 
the present sample from the combination of millimetre and radio 
data. 

Although there might be some fraction of "5.8 jum-peakers" 
that host a certain level of AGN activity, in the majority of these 
sources, the optical to far-lR light is likely dominated by star- 
burst emission. This is confirmed by 1) the high detection rate 
at 1.2mm and the high average value < Fi.2mm >= 1.56 mJy; 
and 2) the correlation between the mm and radio emission 
for the majority of the sources. With star formation rates 
ranging from a few lO^Moyr"' to ~lO^M0yr"' and a mean 
<SFR> > 450-750 Mo yr ', they are powerful ULlRGs wifli 
starburst strengths similar to those of SMGs. 

Significant emission in the mid-lR bands of PA Hs is a 
common featu re in starbu rs ts (see, e.g., Rigopoulou et aL 19991: 
iDesai et al. 20071: iVeiUeux et al. 20091: iFarrah et al. 20091 for 
local galaxies). Although we are still lacking mid-IR spec- 
troscopy for our z ~2 sample to infer the precise strength 
of the objects' PAH emission, there are very good rea- 
sons for thinking that it is strong and comparable to 
that observed in other 24/im-bright starburst galaxies at 
z ~1.5-2 satisfying similar criteria for selection of "4.5 /vm- 
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|3] The black crosses are the sources with S/N < 3 at 1.2 mm. 
The solid black lin e shows the selection criterion for sBzKs of 
iDaddi et al. (2001 . 



" or "5.8^m- peakers" ("Weedman et al. 20061: lYan et al. 20071: 
iMurphv et"ar 2009; Huang et al. 2009' Lonsdale et al. in prep. ) 
or in about 30 "4.5 jum-peakers" (FaiTah et al. 2008). The strong 
PAH emission displayed in these objects accounts for the major 
part of the mid-IR emission in the range ~6-12//m. It is thus 
clear that the PAHs generally contribute much of the flux de- 
tected in the 5/5/fzer/MIPS broad 24 /im band in such "5.8 //m-" 
or "4.5 jum-peakers". As noted, most SMGs have a much weaker 
24fj.m intensity than our sources. It is not yet well understood 
why "5.8/zm-" and "4.5 //m-peakers" happen to have stronger 
PAH emission without a parallel enhancement of the observed 
mm/submm flux density, resulting in a hi gher 24/im/1.2 mm ra- 
tio than the bulk of the SMGs (Lo09). IFarrah et al. (2008 ) sug- 
gest the possibility that star formation is extended on spatial 
scales of l-4kpc in such galaxies. 

In Figure |9l we show the median SED of our sample of 
"5.8jum-peakers", compared to the templates of a starburst, an 
obscured AGN, and an unobscured AGN. It is clear that the me- 
dian observed SED is much better matched by the starburst tem- 
plate than by the obscured AGN SED, where the 5.8 //m bump is 
absent. 

"5.8jum-peakers" (and "4.5 yum-peakers") also have some 
relationship with other broader classes of infrared selected 
high-z star-forming galaxies, as judged by 24 /im inten- 
sity. As already discussed (e.g.. Table 3), there are var- 
ious samples (e.g . iMagUocchetti et al. 20071 : lYan et al. 20051 : 
iHouck et al. 20051: iMurphv et al. 2009h just selected as rel- 
atively strong at 24jum and weak in the optical, espe- 
cially for the purpose of IRS mid-IR spectroscopy. One 
recent example is the so called "Dust-Obscured Galaxies" 
(DOGs) defined as having 5'24um/'^R > 10 00, with variou s 
limits for S24iim such as 300juJy by iDev et al. (20081) : 
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Fig. 11. (Adapted from Fig. 7 of iPaddi et al. 200% . Luminosity 
at 1.4 GHz rest frame vs. luminosity at 8/im rest frame. 
"5.8;um-peakers" from our sample are displayed as large red 
triangles . Data for all 24 yum BzK sources of the GOODS- 
N field (iDaddi et al. 20071) are shown for comparison as small 
black symbols; the three large green squares show the aver- 
age trend vs. 8//m luminosity, including b oth radio-d etected 
and radio-undetected sources from IDaddi et al. (2 007). The 
green stars are the two BzK sources with CO detections from 
IDaddi et al. (20081) . 



Dev & Th e Ndwf s/MIPS Collaboration (20091) . or 100/iJy by 
Pope et a l. (2008a). However, such broad criteria select a mix- 
ture of starbursts and AGN, especially those with power-law 
IRAC SEDs. As seeing from Table 1, practically all sources of 
our sample satisfy or are close to the S2A^lm|SY^ > 1000 flux den- 
sity ratio defining DOGs. 

A better established class of high-z galaxies partially 
selected from op tical-NIR colors is the "BzK" galaxies 
jDaddi et al. 20041) . This includes an "sBzK" sub-class of star- 
forming galaxies and a "pBzK" sub-class of passively evolv- 
ing proto-ellipticals. Fig. [TO] shows that all sources in our sam- 
ple with B, z, and K information are sBzK galaxies. "5.8/im- 
peakers" shar e the same redshif t range as the bulk of BzKs (see 
e.g. Fig. 2 of iDaddi et al. 20()7l) . However, they are at the very 
top of the BzK luminosity function. This is clearly shown in 
Fig. II II where we have put our sam ple of "5.8//m-peakers" on 
top of Fig. 7 of iDaddi et al. (2007h . which displays the rest lu- 
minosity Lv(8/um) (roughly proportional to the bulk of S24fim) 
and LiR (deduced from Li.4ghz) for all BzKs in the GOODS-N 
field. 

We note that the two z ~ 1.5 BzKs where CO was re- 
cently detected by 'Daddi et al. (20081) are "4.5 //m-peakers" with 
SiAfim = 140 and 400 /zJy and Lir close to 10'^ Lq, i.e., compa- 
rable to or below the lower end of the luminosity distribution 
of our sample (Fig. fTTT i. This probably means that CO should 
be easily detectable in most of our "5.8/im-peakers" with the 
current sensitivity of the IRAM Plateau de Bure Interferometer 



(PdBI), as shown by Yan et al. (in prep .), who detected str ong 
CO emission in two "5.8 //m-peakers" of'Sajina et al. (2008"). 

Since sBzK sources are major contributors to star formation 
at z ~2, it is interesting to compare their star formation proper- 
ties with those of our "5.8 /im-peakers". Figure |6] shows that the 
average SFR of "5.8 ^um-pe akers" is signific antly greater than the 
average for sBzKs at z ~2 (IDaddi et al . 2007). Even their SSFRs 
are significantly higher than those of sBzKs with comparable 
masses. This figure confirms that our "5.8yL(m-peakers" represent 
as significant a fraction of the most massive and star-forming 
sBzKs as Fig.[TT]implies. On the other hand, as the average SFR 
of "5.8 /vm-peakers" is significantly lower than those of the bulk 
of classical SMGs, and their stellar mass slightly higher, it is not 
surprising that the mean SSFR of "5.8 //m-peakers" is markedly 
smaller than one of SMGs (Fig.|6]l. 

Another class of powerful starburst galaxi es are the so called 
"submillimetre faint radio galaxies" (SFRGs; iBlain et aL 20041: 
IChapman et al. 2004 12008) . They are defined as radio sources 
with radio fluxes similar to those of SMGs which are not de- 
tected in typical SCUBA surveys. In addition to AGN, this 
class of sources may contain star-forming galaxies character- 
ized by slightly hotter temperatures than typically observed in 
submillimetre galaxies. It is also possible that the radio pro- 
perties and the hotter dust temperatures in these sources might 
be due to the presence of AGN activity, but clear evidence is 
still la cking for the relative fractions of AGN and hotter star- 
bursts ( IChapman et al. 20081: ICasev et al. 2008112009 ). It is thus 
interesting to consider whether our sources undetected at 1 .2 mm 
and lying outside or at the outskirts of the radio- 1.2 mm cor- 
rela tion might belo ng to the class of SFRG. More precisely. 
Cha pman et al. (200 8) propose to reserve the name of SFRG for 
radio sources with Li.4ghz > 10^"* WHz As SFRGs must also 
not be detected in 850 fim surveys with typical sensitivity 4 mJy, 
we assume that SFRGs in our sample may be defined by Li.4ghz 
> lO^^* WHz ' and Fi.2mm < 1.6 mJy. SFRGs are identified with 
these criteria in Table 5. There are 3 confirmed SFRGs and 3 
tentative SFRGs (together with 4 confirmed SMGs and 9 tenta- 
tive SMGs defined as satisfying the limit Fi.2mm > 1.6 mJy, and 
7 sources too weak in radio to be SFRG and at 1.2 mm to be 
SMG). We may thus estimate that among "5.8;um-peakers", at 
least -40% are SMGs, -10-20% may be SFRGs, and -15-20% 
are weak radio and mm sources. 



5.4. Nature of the sources: 2. Weak AGN activity 

Our selection criteria were designed to select starburst 
ULIRGs over AGN through the presence of the redshifted 
1.6 jum stellar "bump". However, similar studies have shown that 
the discrimination against AGN with such criteria is not perfect 
and that a small fraction has a weak but significant AGN con- 
tribution in their mid-IR SED. Since AGNs are powerful NIR 
emitters due to their ability to heat dust up to its sublimation 
temperature, we expect that such emission may smear out the 
1.6//m stellar "bump", flatten i t, or move to long er observed 
wavelengflis dBerta et al. 2007at iDaddi et al. 20071) . However, 
some AGNs can satisfy our selection criteria. Indeed, about 8 
radio-galaxies from the sample of Seymour et al. (2007) follow 
our criteria. But these sources show a sig nificant radio excess 
(I Archibald et al. 200ll:lReuland et al. 2QQ4 due to the AGN. 

To identify, quantify, and characterize any AGN ac- 
tivity which might be present in our sample it would 
be useful to ha ve mid-IR spectroscopic data (see, e.g., 
iFarrah et al. 20081: Weedman et al. 2006), optical/NfIR spec - 
troscopic data (Chapman et al. 2005: iBerta et al. 2007bh . 
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morphological data from high angular resolution radio observa- 
tions (see, e.g., Ric hards et al. 2007; Muxlow et al. 2005; 
Biggs & Ivison 2008), or deep X-ray data (see, e.g., 
Alexander et al. 2005,) . Only limited observations of this 



kind are available for our sample. A spectroscopic observation 
from Kec k/LRIS is available for source L-25 (LH_574364 in 
iBerta et al. 2007b) . Based on one emission line identified as 
Mgll, this source is classified as a type 2 AGN, although its 
radio/1. 2mm flux density ratio is not one of the highest and its 
radio spectral index instead point to a starburst-dominated IR 
SED (TableSi. 

X-ray data of moderate depth (70 ks expo- 
sures) f rom Chandra ar e also available for all our 
sources jPolletta et al. 2006 ^. None of them is detected by 
Chandra to a 0.3-8keV flux limit of 3x10"'^ ergcm"^ s"'. We 
stacked the X-ray images of all 33 sources, but no significant 
detection was obtained. We can only set an upper limit to the 
average 0.3-8 keV flux of 10"'^ erg cm"^ s" 



to an average X-ray luminosity of 2x10^^ ergs Such a low 



which corresponds 

^ergs"' 

signal does not allow us to make any claim about the AGN 
contribution and X-ray properties of our sample. 

The multi-frequency radio observations available for the ma- 
jority of our sources may also allow us to investigate the AGN 
contribution. As discussed in Sec. 4.3, we can find some indica- 
tion of AGN activity based on the 1 .2 mm/radio flux density ratio 
complemented by the spectral index and the radio size; however, 
an extension may as easily reflect well an extended starburst as 
an AGN lobe. 

Two 20 cm radio sources, L-17 and L-11, are specially strong 
(^^20cm >300juJy). The 20 cm flux density of L-17 is slightly 
greater than 300juJy, implying a 1.4 GHz luminosity close to 
10^^ WHz"' (see Fig. [T2l i which is close to the radio loudness 
limit at z ~2 (Jiang et al. 2007; Sajina et al. 2007). L-11, de- 
spite its high 20 cm flux, is not strictly radio loud according to 
this definition. Neither source is resolved in the radio; however. 



this might be not very meaningful, since both sources are in re- 
gions with limited sensitivity at 20 cm. Their spectral indices are 
-0.76 and -0.45, respectively, which is not very discriminating. 
Interestingly, both sources are rather strong and well detected at 
1.2 mm (S/N = 3.5 and 5.3 respectively), implying that, if they 
are AGN, they might experience both starburst and accretion ac- 
tivity. 

The two other sources with evidence of possible AGN-driven 
radio activity based on a high 20 cm/1. 2 mm flux density ratio, 
L-8 and L-14, are not detected at 1.2 mm (S/N < 1 and S/N - 3, 
respectively), L-8 could be an SFRG, and L-14 either an SMG 
or an SFRG. L-8 is also extended in its radio image. 

6. Summary and conclusion 

The aim of this project was to determine the average prop- 
erties of a complete 24 fim flux limited sample of bright Spitzer 
sources selected to be starburst dominated at z ~2, using multi- 
wavelength data. The sample of 33 z ~2 SMG candidates was 
built with all the optically faint sources in a ~0.5 deg^ area of the 
Lockman Hole SWIRE field meeting selection criteria based on 
MIPS/IRAC fluxes. These criteria are F24f2m >400mJy, a peak 
in the 5.8 jum IRAC band due to redshifted 1.6//m stellar emis- 
sion, and r' >23. The J1046H-59 field was selected because 

' Vega 

of the availability of very deep radio observations at 20 cm and 
90 cm with the VLA and at 50 cm with the GMRT. All sources 
in our sample are detected at 50 cm. 

The entire sample has an average 1.2 mm flux density of 
1.56+0. 22 mJy. However, the limited sensitivity allowed only 
four confirmed 4 cr detections, plus nine tentative 3 cr detections. 
Since the average 1 .2 mm flux density, 1 .56 mJy, coiTesponds to 
a 850 jum flux density close to 4 mJy, about half of the sources 
may be considered SMGs. However, their redshifts range from 
z ~ 1.7-2. 3, similarly to the sample in Huang et al. (2009), but 
smaller than the redshift range covered by SMGs. The sample 
selected here is characterized by brighter 24 jum flux densities, 
on average, than those of SMGs, and consequently shows sys- 
tematically lower Fi,2mm/F24^im ratios than classical SMGs. It is 
quite likely that our selection favours the SMGs with the bright- 
est 24 fj.m flux densities, due probably to enhanced PAH emis- 
sion. 

From stacking individual images of the sources, we are able 
to build the median FIR SED of our sample and estimate the 
corresponding Lfir, SFR and Tjust assuming a single tempera- 
ture "greybody" model. The inferred values are T^^st = 37+8 K, 
Lfir =2.5x10'^ Lq, and SFR = 450 Mq yr"' . These estimates in- 
dicate that most of the sources are ULIRGs. However, estimates 
of Lfir for individual sources deduced from the IR-mm SED are 
highly uncertain due to the lack of flux measurements between 
lOOjum and 500 yum. The high quality radio data provide impor- 
tant complementary information on Lfir and the star formation 
rate, since the 1.2 mm/radio flux density ratio of the majority 
of individual sources is consistent with the FlR/radio correla- 
tion, which allows a derivation of Lfir and SFR from the ra- 
dio flux, providing further confirmation that most of the selected 
sources are ULIRGs. The average value of Lfir inferred from the 
FlR-radio correlation is 4. 1 x 10'^ L©; however, this value may be 
overestimated because of an AGN contribution. 

Stellar masses are estimated by modelling the optical-IR 
SED with stellar population synthesis models. They are of or- 
der a few 10" Mq. Roughly scaling with the observed 5.8 //m 
fluxes, they are similar to those of other samples of 24 //m-bright, 
z ~2 Spitzer starbursts, and slightly higher than those of classical 
SMGs. 
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Overall, this sample appears similar to o ther samples of 
Spitzer z ~2 SMGs (Lo09: 'Youn ger et al. 20091) in terms of mil- 
limetre emission, Lfir, and SFR. 

The complete radio detection of all sources provides a good 
estimate of the total star formation rate of such sources. They 
represent a significant fraction of all SMGs in the redshift range 
z ~ 1.7-2.3 (~10-15%). Most of these "5. 8yum-peakers" are star- 
forming BzK galaxies with luminosities at the top of the lumi- 
nosity distribution of sBzKs. 

The surface density of "5.8Aim -peakers" has been found to 
be 61deg"^ bv iFarrah et al. (20061) . This is consistent with 33 
sources in the 0.49 deg^ of our field. We may thus estimate 
that 40-60 similar "5.8//m-peak;ers" per square degree could 
be identified in the full SWIRE survey (49 deg^). Most of them 
should be z ~ 2 starburst ULIRGs. At least half of them may be 
considered to be SMGs, including a small fraction of compos- 
ite obscured AGN/starburst objects. Another significant fraction 
may be considered as SFRGs. 

These results illustrate the power of deep multi-zi studies, es- 
pecially with complete radio data, for analysing populations of 
powerful high-z IR and submm sources. Such deep radio data are 
essential for disentangling starbursts and infrared-bright AGN, 
and for easily providing estimates of their star formation rates. 
We note especially the impressive complete detection of rela- 
tively weak z ~2 SMGs over 0.5 deg^ in a single pointing of 
the GMRT at 6 10 MH z. As already pro ved by the analysis of 
SCUBA sources (e.g., Ilvison et al. 20021) . radio data are essen- 
tial for identifying optical/near-IR counterparts and analysing 
submm surveys. This will be even more crucial for future sur- 
veys at the confusion limits of instruments like Herschel at 300- 
500 yum and SCUBA2 at 850yum. Even as we wait for EVLA 
and the new generation of SKA precursors, our results show that 
the GMRT at 6 10 MHz and even 325 MHz can already currently 
provide sensitivity well matched to wide Herschel surveys. 

It would be interesting to explore further whether the main 
properties which characterize this sample, i.e. strong MIR emis- 
sion, radio activity, and high stellar mass, are related. Some of 
these properties are likely the result of biases introduced by our 
selection; however, this is unlikely to be the case for all of them, 
especially for the radio properties. In particular, a comparison 
between the starburst morphology (traced by young stars, dust, 
PAHs or CO emission, as measured by ALMA or JWST) and the 
radio size would probe whether the radio emission is produced 
by the starburst or by an AGN, and whether the parameters of 
the starburst are different from those of most classical SMGs 
and reveal a different star formation regime. 

We have several multi-wavelength observations planned or 
in progress for this sample to obtain better estimates of redshifts, 
dust temperatures, star formation rates, PAH luminosities, and 
AGN contributions, and thus constrain the dominant emission 
processes, and investigate the evolution and clustering properties 
of these sources. 
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Table 1. Optical, Near-IR and Spitzer mid-IR data of the selected sample 
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1 
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Z4.Z 


OQ ^ 

Zj.j 
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O 1 /I 

Z1.4 
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ZU.Z 


1 n /I 
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3+0.7 


'3/1 Q-Ln 
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y 


QO 0_i_Q 
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J 


Q/i n_i_Q Q 
j4.U± J. J 


/i 1 n o_i_ 1 n 

4iy. /± ly 


J 


T O^ 


C\A7TDCQ Tl 

o W iKni J 1 U4 / 3 


Qo 1 ^n 1 n 1 n 


U 


0/1 1 

Z4.i 


O/I 

Z4.0 


OQ n 

Zj.y 


O'J o 

Zj.Z 


^ OQ /: 

>Zj.O 


on A 
ZU.O 


1 n o 

iy. / 


1 V.U 


44 


4+0.6 


z 1 n_i_n 
J l.U+U 


o 



^O 0_i_Q 

jZ.j± J 


1 

1 


/10 0_i_Q 1 

4/. /± J. J 


000 A_i_ 1 

/ZZ.O+ 1 / 


4 


T OA 


o W iKrlo J 1 U4 /44 
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o 
O 
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>z4.y 


^ O/I n 
>Z4.U 


^ OQ o 

>Zj.Z 


OO z 

Zj.j 


O 1 o 

ZI.J 


on o 
ZU.Z 


ly. j 


36 


3+0.5 


4j.j±U 


o 



/1A 0_i_O 

40.0+Z 


y 


QA n_i_Q Q 
jo.y+j. J 


/I O/I /I _i_ 1 

4y4.4± lo 


c 

J 


L-27 


SWTRF3 1104744 




4 


24.6 


24.5 


23.9 


23.4 


22.6 


21.4 


20.3 


19.6 


27 


8+0.5 


34.6+0 


7 


41.8+3 


4 


<40.0 


522.9+17 


9 


L-28 


SWIRE3J104757 


59+590451 


7 


25.7 


25.9 


25.1 


24.4 


23.2 


21.1 


20.3 


19.3 


32 


6+0.6 


40.7±0 


8 


48.5+3 


2 


<40.0 


418.7+18 


8 


L-29 


SWIRE3J 104807 


86+591119 


5 


24.6 


24.8 


24.2 


23.5 


22.5 


20.8 


19.8 


19.1 


40 


7+0.6 


49.3+0 


8 


64.2+3 


3 


38.6±3.2 


407.2+19 


1 


L-30 


SWIRE3J104816 


81+584658 


5 


24.1 


24,7 


24.3 


23.6 


>23.6 








35 


6+0.6 


45.7+0 


9 


49.1+3 


3 


<40.0 


469.9+16 


4 


L-31 


SWIRE3J104830 


71+585659 


3 


>24.5 


>24.9 


>24.0 


>23.2 




>24.2 


22.8 


21.0 


17 


4+0.5 


24.2+0 


9 


43.8+3 


5 


<40.0 


639.9+15 


6 


L-32 


SWIRE3J104841 


00+590056 


8 


>24.5 


>24.9 


>24.0 


>23.2 




22.9 


20.9 


20.0 


27 


0+0.6 


37.1+0 


7 


44.2+3 


3 


30.6+3.0 


406.2+19 


3 


L-33 


SWIRE3J104848 


23+585059 


3 


20.1 


25.1 


24.3 


23.8 


22.9 


21.4 


20.3 


19.7 


32 


1+0.7 


37.6+0 


9 


51.4+3 


7 


33.4+3.8 


559.7+20 


3 



The optical magnitude limits correspond to 90% completeness values derived in the field where our sample was selected, but note that the coverage of the optical observations is not uniform 
the field. 

The optical data come from NOAO and the Near-IR magnitudes from UKIRT (see text) 
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Table 2. Related Spitzer-sslsctsd z ~ 2 ULIRGs samples. 



Sample 


N 




< ^24 //m > 


IR color selection 


r 


density 












Vega 


deg"^ 


This work 


33 


>400 


566 


-f3.6;im<-f4.5/im<-'^5.8;jm>-f8.0//m 


>23 


65 


Lonsdale et al. 2008 


61 


>400 


819 


-f3.6//m<f4.5//ra<-f5.8;jm>f"8.0(/ra 


>23 


55 


Farrah et al. 2008 


32 


>500 


726 


^3.6;im<-f'4.5//ra>-f"5.8;jm and 5^n,>Fg fl^^ 


>23 


34 


Huang et al. 2009" 


12 


>500 


689 


0<[3.6]-[4.5]<0.4 
-0.7<[3.6]-[8.0]<0.5 






Magliocchetti et al. 2007 


793 


>350 






>25.5 


200 


Yan et al. 2005 


52 


>900 




log,o(y/,(24y^m)/y/,(8A/m))>0.5 







logio(y/,(24//m)/y/,(0.7Aim))> 1 .0 



" The sample of lHuang et al. (2009h differs by only one source from the one o f lYounger et al. (20091) . 



Table 3. Derived values from the stacked flux densities". 



Sample 


N 


z 


Fl.lmm 




^1 .2 mm/ ^24 fiui 


L]A GHz 


/ 

'-'FIR.mm 


^FIR.radio 


SFRmm 


SFRradio 


q 








(mJy) 


0"Jy) 




(lO^-'WHz-') 


(10'2 




(Mo 


yr-1) 




All 


33 


2.08 


1.56+0.22 


92.6±12.0 


2.76+0.50 


2.03+0.26 


2.45+0.35 


4.14+0.53 


441 


745 


2.11+0.12 


S/N>3* 


13 


2.11 


2.83+0.14 


128.5+26.0 


4.94+0.51 


2.89+0.58 


3.99+0.20 


5.90+1.18 


718 


1062 


2.17+0.11 


S/N<3* 


20 


2.06 


0.72±0.18 


69.3±7.4 


1.28+0.40 


1.48+0.16 


1.38+0.35 


3.02+0.32 


248 


543 


2.00+0.16 



"The average stacked flux densities are computed with equal weights to avoid biases. Their approximate rms are computed as the standard 
deviation of the mean. 
'S/N at 1.2 mm. 

'^ipiR.mm is computed with (Eq. 3) and the dust temperatures coming from the stacked SED (see text): 37 K for the entire sample, 36 K for the 13 
sources with S/N > 3 at 1.2 mm, and 39 K for the 20 sources with S/N < 3. 
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Table 4. 1 .2 mm and radio data. 



Source ID 


f 1.2mm° 

(mJy) 


^"20 cm" 


FsOcm" 

(MJy) 


Fgocm" major axis 

O^Jy) (") 


minor axis 
(") 


a" 


q 


1.2 mm detections S/N>4 


L-1 


2.95+0.66 


76.3+9.0 


142+19 


327±74 


<2.9 


X 


-0.84±0.18 


9 40+0.20 
^•^^-0.20 


L-9 


4.00±0.55 


116.5±9.2 


195 ±22 


266±69 


1.9 


0.7 


-0.54±0.03 


7 a c +0.20 
■-^-^-0.20 


L-11 


3.08+0.58 


314.8±4.1 


425±43 


651±72 


X 


X 


-0.45±0.06 




L-22 


3.41±0.73 


100.9±7.1 


238±26 


447±70 


<1.7 


X 


-1.02±0.01 




1.2 mm tentative detections S/N>3 


L-5 


2.75±0.76 


71.0±20.0 


177±23 


322±74 


<4.5 


X 


-1.03±0.06 


2.28-;- 


L-14 


2.13±0.71 


159.5±6.0 


322±33 


426±70 


0.7 


X 


-0.73±0.09 


1 1 +0.26 


L-15 


2.36+0.62 


68.6±3.7 


181±21 


258±69 


<1.0 


X 


-1.03±0.14 


9 Ti +0.20 
^■-'-'--0.20 


L-17 


2.24+0.64 


341.0+33.0 


756±78 


1052±75 


X 


X 


-0.76±0.10 


1 ftQ+0.20 
i-°°-0.20 


L-20 


2.66±0.78 


51.1±4.7 


166±19 


264±72 


<1.5 


X 


-1.25±0.18 


1 95+"" 


L-21 


3.09±0.81 


138.3±35.8 


176±38 


350±76 


4.1 


1.8 


-0.64±0.28 


2 51+0.37 


L-23 


3.13±0.86 


86.4±7.2 


172±21 


242±71 


<1.3 


X 


-0.77±0.07 


9 7 Q +0.23 
• -^-0.23 


L-25 


2.56±0.74 


69.0±9.0 


166±20 


348±72 


<1.4 


X 


-1.09±0.04 


2 63+"-^5 


L-27 


2.48±0.73 


77.2±13.7 


106±18 


<219 


<2.1 


X 


-0.38±0.39 


9 47+0.27 
^•^'-0.27 


Other non detections at 1.2 mm 


L-2 


<3.04(1.46±0.79) 


<42.0 


90±17 


<215 


X 


X 


[-0.64] 




L-3 


<3.33 (1.04±1.11) 


105.0±35.0 


50±17 


<227 


<7.6 


X 


+0.89±0.85 




L-4 


<3.42 (0.76+1.14) 


<93.0 


76±17 


<220 


X 


X 


[-0.64] 




L-6 


<2.22 (0.59±0.74) 


37.1±6.4 


69±12 


<202 


<1.4 


X 


-0.75±0.42 




L-7 


<2.97 (1.39±0.79) 


101.7±13.3 


200±23 


<216 


2.0 


0.8 


-0.81±0.29 




L-8 


<2.70 (-0.36±0.90) 


124.0±14.0 


219±25 


250±72 


1.4 


X 


-0.58±0.13 




L-10 


<3.01 (1.39±0.81) 


67.2±11.7 


94±17 




<3.7 


X 


-0.40±0.37 




L-12 


<3.64 (2.06±0.79) 


94.4+4.7 


182±28 


315±69 


<1.3 


X 


-0.81±0.02 




L-13 


<3.15 (-0.49±1.05) 


<39.0 


58±16 


<219 


X 


X 


[-0.64] 




L-16 


<2.58 (0.51 ±0.86) 




145±32 


319±77 


X 


X 


-1.25±0.76 




L-18 


<2.43 (-0.72±0.81) 


93.2±9.4 


162±19 


<209 


2.4 


X 


-0.67±0.36 




L-19 


<2.91 (0.89±0.97) 


102.2±12.0 


138±18 


231±69 


3.9 


1.8 


-0.44±0.13 




L-24 


<2.64(0.14±0.88) 


<24.0 


45±14 


<211 


X 


X 


[-0.64] 




L-26 


<2.58 (0.37±0.86) 


54.9±13.4 


84±14 


255±72 


3.5 


X 


-0.86±0.53 




L-28 


<2.96(1.64±0.66) 


<23.0 


47±12 


<213 


X 


X 


[-0.64] 




L-29 


<2.94 (0.45±0.98) 


55.0±12.0 


76±16 


<216 


X 


X 


-0.39±0.52 




L-30 


<1.89 (0.42±0.63) 


<70.0 


84±18 


<236 


X 


X 


[-0.64] 




L-31 


<2.82 (1.62±0.60) 


121.0±39.0 


157±20 


<221 


6.6 


X 


-0.31±0.50 




L-32 


<2.89 (1.61±0.64) 


<48.0 


60±16 




X 


X 


[-0.64] 




L-33 


<3.27 (-0.28+1.09) 


109.0±33.0 


87±19 


<235 


<10.0 


X 


+0.27±0.61 





"The limits on Fi 2mm are '3cr' limits, except for the sources with 1.5 < S/N < 3 where they are the observed value plus 2 cr. The values in brackets 
are the observed values used in stacks. 
*The radio flux density limits are '3cr' limits. 

■^The values of a in brackets show the cases where a cannot be determined and is assumed equal to -0.64 (see text). 
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Table 5. Luminosities and star formation rates. 



Source 


Zphol 


I 1' 

1.4 GHz 


^FIR.mm 


^FlR,radio 




vL,,(1.6yum) 


M* 


SMG/SFRG 






(lO^^'WHz-') 


(10'- 




(Mo yr- 




(10" Lo) 


(10" Mo) 




L-1 


1.92 


1.60 


4.63+1.04 


3.27+0.81 


834 


588 


1.67 


1.03 


SMG 


L-2 


1.92 


0.96 


<4.77 


1.96+0.39 


<859 


353 


1.61 


1.08 


weak Radio? 


L-3 


1.73 


0.33 


<5.22 


0.66+0.52 


<941 


<284 


1.66 


1.13 


weak Radio 


L-4 


2.13 


1.02 


<5.37 


2.08+3.08 


<966 


<1664 


3.02 


2.85 


SFRG? 


L-5 


2.04 


2.40 


4.32±1.19 


4.90+0.92 


111 


883 


2.53 


1.36 


SMG? 


L-6 


2.26 


1.10 


<3.49 


2.25+0.79 


<627 


405 


2.07 


1.07 


SFRG? 


L-7 


1.89 


2.15 


<4.66 


4.35+1.06 


<839 


788 


2.19 


1.55 


7 


L-8 


2.42 


3.80 


<4.24 


[7.75+1.09] 


<763 


[1396] 


2.89 


[1.52] 


SFRG 


L-9 


2.33 


2.99 


6.28+0.86 


6.11+1.94 


1130 


1100 


2.71 


2.66 


SMG 


L-10 


2.45 


1.56 


<4.73 


3.19+1.05 


<851 


574 


4.01 


2.31 


7 


L-11 


1.92 


4.52 


4.84+0.91 


[9.21+3.63] 


870 


[1658] 


1.57 


[0.78] 


SMG 


L-12 


1.71 


1.56 


<5.71 


3.19+1.12 


<1029 


575 


3.19 


3.09 


7 


L-13 


2.22 


0.86 


<4.95 


1.75+0.89 


<891 


315 


1.60 


0.87 


weak Radio 


L-14 


2.36 


5.53 


3.34+1.11 


[11.30+2.99] 


602 


[2029] 


3.37 


[1.74] 


SMG? 


L-15 


1.86 


1.92 


3.71+0.97 


3.93+0.76 


667 


707 


3.09 


3.21 


SMG? 


L-16 


1.73 


1.36 


<4.05 


2.77+1.51 


<729 


499 


1.34 


0.92 


SFRG? 


L-17 


1.90 


7.90 


3.52+1.00 


[16.12+2.70] 


633 


[2902] 


2.27 


[1.07] 


SMG? 


L-18 


1.85 


1.60 


<3.82 


3.27+0.86 


<688 


589 


3.05 


2.17 


SFRG 


L-19 


1.69 


1.08 


<4.57 


2.21+0.53 


<822 


398 


1.67 


0.99 


7 


L-20 


2.69 


4.98 


4.18+1.22 


10.20+0.95 


752 


1830 


5.68 


7.16 


SMG? 


L-21 


1.78 


1.61 


4.85+1.27 


3.23+1.94 


873 


592 


1.74 


1.12 


SMG? 


L-22 


2.57 


5.70 


5.35+1.15 


11.62+4.323 


964 


2092 


4.11 


3.06 


SMG 


L-23 


2.38 


3.08 


4.91 + 1.35 


6.28+1.63 


885 


1130 


2.27 


1.27 


SMG? 


L-24 


2.07 


0.57 


<4.14 


1.16+0.29 


<746 


<209 


1.82 


1.02 


weak Radio 


L-25 


1.47° 


1.01 


4.02+1.16 


2.07+0.59 


723 


372 


1.29 


0.84 


SMG? 


L-26 


1.72 


0.73 


<4.05 


1.49+1.07 


<729 


<579 


1.55 


1.57 


weak Radio 


L-27 


2.20 


1.41 


3.89+1.15 


2.87+0.96 


701 


517 


2.18 


1.21 


SMG? 


L-28 


1.88 


0.48 


<4.65 


0.98+0.20 


<836 


<176 


1.79 


1.29 


7 


L-29 


1.93 


0.77 


<4.62 


1.58+0.68 


<831 


284 


2.33 


1.62 


weak Radio 


L-30 


2.38 


1.45 


<2.97 


2.96+2.39 


<534 


<1290 


3.17 


2.01 


SFRG 


L-31 


2.74 


3.13 


<4.43 


6.39+3.43 


<797 


1150 


3.65 


1.37 


7 


L-32 


2.38 


1.03 


<4.53+1.00 


2.10+1.54 


<815 


<278 


2.77 


1.84 


7 


L-33 


2.18 


0.92 


<5.13 


1.87+0.97 


<924 


337 


2.47 


1.49 


weak Radio? 



"Spectroscopic redshift from lBerta et al. (2007bl) 

'The value of Li.4ghz! the rest-frame luminosity at 1.4 GHz, is the weighted average of the two determinations of Li 4ghz from 20 cm and 50 cm 
flux densities when both are available. 

^The values of LpiR.radio, the infrared SFR derived from the radio flux densities and the stellar mass are dubious for sources with high 

20 cm/1.2 mm flux densities ratios. They are reported, in brackets. 

''SFR calculated by Kennicutt's formula (Eq. 6) with a Salpeter IMF (see text). 

The upper limits of LpiR,,™,, and SFR,n,^ are '3cr' limits, except for the sources with 1.5 < S/N < 3 where they are the observed value plus 2 a. 
The upper limits of LpiRjadio and SFRradio are '3cr' limits. 
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Table 6. Results from the co-added MIPS images. 



Sample 


N 








F\bOfiml< ^1.2miTi >" 










from stack median (mjy) (mjy) 






GO 30391 data 


All 


21 


571±27 


2.78+0.33 


13.47+2.76 1.67+0.27 


8.1+3.0 


4.9+0.8 


1.2mm SNR>3 


10 


585+32 


2.56+0.38 


14.63+4.05 2.77+0.18 


5.3+1.8 


4.4+0.9 


1.2 mm SNR<3 


11 


559+44 


3.08+0.67 


11.95+4.08 0.68+0.24 


17.6+12.2 


5.5+1.6 


Younger et al. (2009) data 


All 


12 


680+50 


3.3+0.6 


21.7+7.9 1.4+0.3 


15.5+9.0 


4.9+1.3 



°Fi60fim/< Fi2mm > and F^a^J < F24f,m > are computed with the stack medians for Fjo^ and Figo^, 
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Appendix A: SED fits and photometric redshifts 

The optical (U^ri z), NIR (JHK), and MIR (3.6-24jum) SEDs of each source have bee n fitted with a library of 18 star-forming 
galaxy templates ("Poll etta et al. 2 007*) using the Hyper-z code ("Bo lzonella et al. 20001) and the same procedure described 
in [Lonsdale et al. (200C 

In Fig. lA.ll we show the optical-IR SED and the best-fit template and corresponding photometric 
redshift of each source. In seven cases we also report a secondary solution corresponding to a second minimum in the;^'^ 
distribution if associated with a different template than the primary solution. In case a spectroscopic redshift is available (source 
L-25), the best-fit template at the spectroscopic redshift is also reported. 
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Fig. A.l. Optical-IR SEDs of our sources. The solid curves represent the best-fit template of the optical-MIR data (up to 24 yum). The 
dashed blue curves correspond to fits with increasing The template name and photometric redshifts of each fit are annotated. For 
L-25, the dotted curve corresponds to the best-fit template plotted at the spectroscopic redshift. The source ID number is reported 
in each panel. Downward arrows represent 5cr upper limits at optical and infrared wavelengths. 
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g.A.l. Continued 
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